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We evaluate the effective number of atoms in experiments wliere a 
probe laser beam with a Gaussian transverse profile passes through an 
atomic medium consisting in a cold atom cloud released from a magneto- 
optical trap. Considering the case where the initial distribution is a Gaus- 
^ ' sian function of position and of velocity, we give a quantitative description 

, of the temporal variation of the effective number while the cloud is ex- 

1/^ . ploding and falling down. We discuss the two cases where the effective 

' number is defined from the linear and nonlinear phaseshifts respectively. 

We also evaluate the fiuctuations of the effective atomic number by cal- 
■ culating their correlation functions and the associated noise spectra. We 

finally estimate the effect of these fiuctuations on experiments where the 
probe beam passes through a cavity containing the atomic cloud. 
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_ 1 Introduction 

, In recent years, the development of magneto-optical traps |l|] has opened the 

way to the use of cold atoms for various purposes such as atomic interferometry, 
and nonlinear or quantum optics. In many of these experiments, a cloud of cold 
atoms is first trapped and then released when the laser beams and the magnetic 
field of the trap are turned off. The experiment is thus performed while the cloud 
}J] ' is undergoing a ballistic expansion due to the velocity distribution in the trap 

5^ \ and falling down under the influence of gravity. Models using a Monte-Carlo 

simulation of the evolution of the atomic positions and velocities have been 
developed to analyze the free fall of the cloud of cold atoms. They have been 
used to obtain informations about the initial position and velocity distribution in 
the trap from fluorescence measurements, particularly to evaluate temperature 
in the trap when it is released and to investigate deviations from a Maxwellian 
velocity distribution (||] and references therein). 

Here, we want to concentrate on the category of experiments where the 
falling cloud of cold atoms is used in spectroscopy, nonlinear or quantum optics 
experiments. In this case, the cloud interacts with a probe laser beam that 
passes through it. A crucial point for these experiments lies in the definition 



1 



of an effective number of atoms contributing to the interaction with the probe 
beam and also to variation of this effective atomic number as a function of time. 

A first aim of the present paper is to provide a definition of the effective 
number of atoms in the probe beam, as well as a precise modelisation of its 
temporal evolution. We will show that this is a quite involved problem in the 
general case, and in particular that different definitions may be given, depending 
upon the specific physical effect chosen as a reference for the evaluation of the 
atomic number. We will discuss the example where the linear phase produced 
by the atomic medium is chosen as a reference, as well as the case where the 
nonlinear phase shift is considered. This general discussion will lead to simpler 
conclusions in the two limiting cases where either the transverse size of the probe 
beam is much smaller than the cloud size or the Rayleigh length characteristic 
of the divergence of the probe beam is much larger than the cloud size. 

A further problem is that this effective number of atoms has fluctuations. 
These fluctuations must be added to the fluctuations usually evaluated with the 
number of atoms kept constant and they may have an influence in spectroscopy, 
metrology or quantum optics experiments. In this class of experiments indeed, 
the signal is often measured with such a high precision that any excess noise 
produced by the atomic medium can appear as a limiting factor for the sensi- 
tivity. 

This is why in a second part of this paper we will evaluate atomic number 
fluctuations by calculating their characteristic correlation functions. We will 
show in particular that it is possible to calculate time dependent noise spectra, 
when the transverse size of the probe beam is much smaller than the cloud 
size. In this case indeed, the correlation time, expected to be of the order of 
the time of flight through the probe beam, is much smaller than the typical 
variation time of the mean number of atoms, of the order of the time of flight 
through the cloud, so that the effective atomic number may be considered as a 
quasistationary random variable. 

The theoretical developments of the present paper are in particular essential 
for the understanding and interpretation of experimental results we obtained 
recently The experiment consists in the measurement of quantum fluctua- 
tions of a weak probe beam passing through a cloud of cold atoms expanding 
and falling down after it has been released from a standard magneto-optical 
trap. To enhance the interaction between atoms and light, the atoms are placed 
inside an optical cavity. The variation of the number of atoms leads to a con- 
tinuous variation of the linear phaseshift of the laser beam, which is sufhcient 
to scan the cavity across resonance. Due to the non-stationary character of the 
measurement of quantum fluctuations, the knowledge of the time variation of 
the number of atoms plays a crucial role in the interpretation of experimental 
results. The fit with the experimental results, using the theory developed in the 
present paper, will be presented in a further publication While a variation 
of the mean atomic number scans the cavity detuning through a variation of 
the linear phase shift, atomic number fluctuations result in fluctuations of the 
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detuning which have an effect on the fluctuations of the probe beam which are 
measured in the experiment. If they are large enough, they may enter exponen- 
tially in the expression for the field fluctuations after interaction, and therefore 
produce excess noise not only at the frequencies contained in their own noise 
spectrum, but also at higher frequencies through a multiplicative noise process- 
ing 1^. The evaluation of atomic number fluctuations contained in the present 
paper are particularly useful to delineate under which conditions these fluctua- 
tions degrade squeezing measurements. 

As already pointed out by Weiss et al. 1^, the density distribution of a 
freely falling atom cloud can be calculated analytically at any instant, if a 
Gaussian profile in the trap is assumed for position density as well as for velocity 
distribution. Experiments by Drewsen et al. ||] have shown this assumption to 
be justified if multiple scattering effects become negligcable, that means for 
large clouds and atoms with a temperature in the 100 /iK to mK range 0. 
These conditions are met in the experiment that we are modelizing jsj. When 
the number of trapped atoms is deduced from their fluorescence signal, the 
detectivity function does not have a Gaussian profile and does in general not 
lead to analytical results. In the case studied in this paper, the experimental 
determination of the number of atoms relies on a phase shift measurement of 
a probe laser beam. It may thus be expected that the detectivity function is 
given by the Gaussian profile of the laser beam, leading to analytical solutions 
for the mean atomic number and for correlation functions of atomic number 
fluctuations. 

In Section 2 of this paper, we describe our model in more detail and show 
how the effective number of atoms is deduced from phase shift measurements, 
accounting for the transverse mode structure of the probe beam. From this, we 
analyze in Section 3 the expression of the case of a linear phase shift measure- 
ment, considering the general case of an atomic motion under gravity. In Section 
4, we address the problem of a saturated phase shift. The correlation functions 
and the noise spectrum of atomic number fluctuations are then calculated in 
Section 5. When evaluating fluctuations, we restrict our attention to the partic- 
ular case of the atomic number deduced from the linear phase shift since we are 
mainly interested in potential effects of linear phase fluctuations. We finally give 
an estimate of the effect of the latter on the probe-cavity detuning in Section 6. 

2 Description of the model 

As already discussed, we consider that the number of atoms is monitored through 
a phase measurement performed on a probe laser beam. Hence, this number will 
be obtained as an integral of the atomic density over the volume, with a weight 
function fitting the Gaussian profile of the laser beam. 

To determine precisely this weight function, we introduce the Gaussian mode 
u (r) (r stands for the three-dimensional position (x, y, z)) describing the intra- 
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cavity field propagating along the x direction and coming onto the atoms; 



u{t) 



1 



■ exp 



IT W (x) 

w (x) is the position-dependent beam size: 

{x) = Wq [1 



Ir 



w"^ {x) 



iip (r) 



A 



(1) 



(2) 



(3) 



where u>o is the beam waist, If; the Raylcigh divergence length and A the laser 
wavelength. The field phase (p (r) is the sum of contributions representing prop- 
agation phase and curvature of the wave at longitudinal position x: 



ip{r) = — h arctan 

A 



TT {y^ + Z^) X 
A + 



The squared modulus of the Gaussian mode: 

/(r) 



\uir)\^ = 



S{x) 



(4) 



(5) 



may be written in terms of a Gaussian weight function / (r) normalized to unity 
on the beam axis and an effective beam section S: 



f (r) = exp 



Six) 



w"^ (x) 

'KVp' {x) 



(6) 
(7) 



The normalisation of the Gaussian mode u (r) has been chosen so that: 

\u{r)\^dy dz = l (8) 



This relation holds for any position x along the propagation axis. 

We will write the incident field seen at time t by the atoms located at point 
r in the interaction representation with respect to the laser frequency ujl'- 



E {r,t) = e-^'^^E (r,i) + e^'^^E (r,t)^ 



(9) 



where E {r,t) is the product of the Gaussian profile u (r) by a mode amplitude 
Ait): 



Eiv,t) 



Tiu>i 
2eoc 



Ait)uiv) 



(10) 
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This definition is such that A{t)'' A {t) is the number of photons passing through 
a beam section per unit time. 

As explained below, the field radiated by the atoms will have a different 
transverse profile as a consequence of the nonlinear character of light-matter 
interaction. The Gaussian mode u (r) is the lowest order mode mqo of the family 
of orthogonal Hermite modes. We will consider that the higher-order Hermite 
modes Umn (with m ^ or n 0) are not resonant at the same frequency as 
the lowest order mode Uqq. Hence, any field radiated by the atoms into these 
higher-order modes will not be efficiently coupled back onto the atoms and can 
therefore be disregarded. We will then calculate the field amplitude radiated 
into the fundamental mode by a mere projection. 

The effect of an optically thin atomic layer of length dx may be described 
by a local modification dE (r)of the field: 

E{r)^ E (r) + dE (r) (11) 

with §: 

Q\2 

dE{r)^-—p{r)a{r)E{r)dx (12) 
An 

Note that we have treated the motion of atoms as quasistatic and neglected any 
consequence of the Doppler effect. As a result, we have used the atomic density 
p (r) integrated over the velocity distribution, a (r) is the atomic polarizability 
measured as a dimensionless number: 

The linear polarizability ai is a function of the dimensionless detuning d between 
laser frequency lol and atomic resonance frequency loq normalized to the decay 
rate 7 of the atomic dipole: 

= ihs ^''^ 

^^u^o^ (15) 

7 

The saturation parameter s is proportional to the local laser intensity |/3 (r)|^ 
measured as a dimensionless number: 



doE{v) /3A2 A 
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where do is the matrix element of the atomic dipole. 

As a first step, we consider that the saturation of polarizabifity may be 
neglected; we will come back to the nonlinear case later on. The number of 
atoms is then evaluated by monitoring the linear phase shift produced by the 
atomic cloud, with dE (r) written as: 

3\2 

dE (r) = — —p (r) aiE (r) dx (18) 
47r 

Even in this linear case, the radiated field dE (r) docs not have the same spatial 
variation as the incident field E (r) when the atomic density p (r) is position- 
dependent. Ignoring the field radiated into higher-order modes, one gets the 
modification of the field amplitude dA in the Gaussian mode by projecting 
dE (r) onto u (r): 

^ —dA = [ dE (r) u* (r) dy dz (19) 

Assuming that the modification of the field by the whole atomic cloud remains 
small (assumption of an optically thin atomic cloud), we obtain this modification 
by summing up the contributions of all the atomic layers as: 

- = --a.a (20) 

with : 

a = J\u{r)f p{r) dr (21) 
The modification of the field intensity is: 

d (A^A) =AUA+ (dA^ A = -A^A^ — ^ (22) 

27r 1 + 

As the absorption cross section is (3A^/27r) at resonance and is divided by 
(l + (5^) otherwise, this equation shows that a must be interpreted as the ef- 
fective number of atoms present in the detection beam per unit beam section. 
This interpretation may be made more precise by rewriting a as: 



dn (x) 



(23) 



where dn (x) is the number of atoms present in the detection beam in a thin 
layer of length dx, that is the number obtained by integrating the atomic density 
p (r) over the beam area with a Gaussian weight function / (r): 

^= j f{v)p{v)dydz (24) 
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This effective number has to be distinguished from the number of atoms present 
in the whole cloud in a thin layer of length dx, which is obtained by integrating 
the atomic density p (r) : 



dN{x) 



p (r) dy dz (25) 



dx 

The atomic medium considered in this paper consists in a freely falling cloud 
of cold atoms. We will assume that the phase space distribution 7r(r, v,0), at 
time i = chosen as the beginning of the free fall, is a Gaussian function of 
position as well as of velocity: 

Af is the total number of atoms while measures the radius of the initial cloud 
and (7y the thermal velocity related to the trap temperature T through: 

= ^ (27) 

m 

As already mentioned in the introduction, this Gaussian approximation gives a 
rather good description of the trap distribution [|[ It presents the great ad- 
vantage that it will lead to analytical expressions for the integrals encountered 
in the evaluation of the mean atomic number and of atomic number fluctua- 
tions, at least in the case where the Rayleigh divergence length is larger than 
the cloud size. For the sake of simplicity, the initial cloud has been considered to 
be centrered on the detection beam, and the detection beam supposed to prop- 
agate in horizontal direction. Note however that releasing these two simplifying 
assumptions would not change the Gaussian character of the atomic density, 
which is the key point of the following calculations. 

The motion of atoms at later times t is determined by the laws of free fall, 
which preserve the phase space volume: 

dr (t) dv (t) = dr (0) dv (0) (28) 

The time-evolution of the phase space distribution tt (r, v,i) is deduced from the 
conservation of the number of atoms (Liouville's theorem): 

n(^r + vt+ igt^, V + gt, ?j ^TT (r, v,0) (29) 
with g oriented downwards in the z direction, or equivalently: 

TT (r, v,t) = TT (^r - v< + igt^, V - gt, 0^ (30) 

The atomic density is defined as the integral over velocity of the phase space 
distribution: 
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J vr (r, v,t) dv 



(31) 



This Gaussian integral is easily evaluated 



3/2 



exp 



2(a2+a2i2) 



(32) 



For example, the atomic density at the initial center of the cloud (r = 0) is: 



[27rt72]3/2 



t2 



3/2 



exp 



-2 (r2 + f2) 



with: 



2y2cr„ 



(33) 



(34) 



(35) 



This shows that the atomic cloud undergoes a ballistic expansion on a time scale 
Tr, the time of flight through the trap radius Ur for an atom flying at mean 
thermal velocity CTu, while at the same time the cloud is falling down under the 
influence of gravity on a time scale Tg. When the time scale Tg is longer than 
Tr, the effect of gravity may be disregarded for short times {t <^ Tg) and the 
atomic density varies mainly as a consequence of ballistic expansion. Gravity 
leads to an exponential decrease at longer times {t ^ Tg) and the variation of 
the atomic density may be rewritten in this case: 



p(0,i) 



[2t. 



i3/2 



i2 



3/2 



exp 



t 



2 1 



(36) 



A third time scale will be found in the following, which will measure the time 
of flight through the detection beam: 



w {x) 
2a r, 



(37) 



3 Effective number of atoms deduced from the 
linear phase shift 

Using the results of the previous section, the linear phase shift of the field due 
to the interaction with the atomic cloud may be written in terms of equations 
20i I23, P4) where we have introduced a time-dependence of the effective atomic 
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number per beam section a (t) in order to account for the quasistatic time- 
variation of the atomic density p (r, t) : 



_^.a,—ait) (38) 
.(.)^/^ (39) 



dn (x, t) 

dx 



= J !{v)p{v,t)dydz (40) 



The latter integral over the transverse variables y and z has a Gaussian character 
and can be readily performed. The resulting expression for the number of atoms 
dn{x, t) present at time t in the detection beam in a thin layer of length dx may 
be written: 



dn (x, t) 



dU{x,t) w2 [x) 



4 ((72+^2*2) 



■ w 



Hx] 



exp 



(41) 



dM {x, t) is the number of atoms present at time t in the whole cloud in a thin 
layer of length dx (cf . definition ^5|) : 



dN {x,t) 
dx 



V'27r(a2 + a2t2) 



exp 



2(a2+a2t2) 



(42) 



The integral over the longitudinal variable x keeps a Gaussian character either 
when the beam size is much smaller than the trap radius (w <C CTr) or when the 
Rayleigh divergence length is much larger than the cloud size ((Tr ^ Ib)- We 
discuss these two limiting cases in the next paragraphs, and recall that both 
assumptions are fulfilled in the experimental situation of reference 

We consider first the limiting case w <^ Ur, but with OrjlR arbitrary, where 
the number of atoms dn{x,t) given by equation ( pl| ) is proportional to {x), 
and therefore also to the beam section S{x). The x-dependent factor S{x) dis- 
appears in the expression (|3^) of a (t) which is thus a Gaussian integral. Its 
evaluation leads to: 



a{t) 



2^2 (T2+t2) 



exp 



(r2+t2) 



(43) 



where and Tg are the already introduced time scales which correspond re- 
spectively to free flight through the trap and free fall (compare with equation 

In the case where the Rayleigh divergence length is much larger than the 
cloud size (ct^ <SC Ir with w/ur arbitrary), the x-dependence of the beam sec- 
tion S can be ignored and a (t) is again a Gaussian integral. For an arbitrary 
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beam size w, the time of flight through the probe beam becomes a relevant 
parameter. After the evaluation of the integral, we obtain a {t) as: 



a it) 



exp 



(44) 



where is the time of flight through the probe beam. When both assumptions, 
a large Rayleigh divergence length and a small beam size {w <ti Ur and cr^ ^r), 
are valid, expression (Q) is identical to expression (^3|) which is valid for any 
ratio Ot/Ir- 

As already discussed with the help of equation (p3|), the effect of gravity 
may be disregarded for times t <ti Tg when Tj. <C Tg. It follows that n(t) thus 
varies as a Lorentzian function of time, as a consequence of ballistic explosion. 
Gravity leads to an exponential decrease at longer times. In the particular case 
Tu, <C Tj. <C Tg for example, these various regimes of time-dependence of a (t) 
may be described in a single formula: 



a{t) 



exp 



t2 



(45) 



This condition <§; Tg corresponds to a high temperature limit, since it may 



be rewritten from equations (2734 35) 



kBT-» 



mar |g| 
2^/2 



(46) 



In the opposite low temperature limit, the variation of cr [t) is mainly determined 
for all times by the exponential factor which describes the free fall of the cloud 
under the effect of gravity. 



4 Effective number of atoms deduced from the 
nonlinear phase shift 

In the previous section, we have studied the field modification in a linear regime, 
and deduced an effective number of atoms from a linear phase shift. We now 
come to the general case where the atomic polarizability may be saturated. 

The transverse variation of dE (r) is then completely different from the trans- 
verse variation of the incident field E (r) as a consequence of the dependence of 
the saturation parameter (see equations |l2[]T^): 



s (r) = s„ {x) f (r) (48) 
Sm is the saturation parameter evaluated on the beam axis: 
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As in the linear case, we ignore the field radiated into higher-order modes and we 
obtain the modification of the field amplitude dA projected onto the Gaussian 
mode through equation (|l9|): 

^ = (50) 

where Us is now modified by saturation: 

-W^/^ (51) 

In order to evaluate Us , we will first expand dug in powers of the saturation 
parameter: 

oo 

dris {x, t) = Y^ (-2s„ {x)f dn^^+^'> {x, t) (53) 



fc=0 



where: 



= J f{ryp{v,t)dydz (54) 



dx 

dn^^''{x, t) is exactly the expression (^0|) of the effective number of atoms dn{x^ t) 
in the detection beam in a thin layer of length dx, as it was evaluated in the 
previous section. Due to the Gaussian shape of function /(r), function /(r)-' has 
the same expression as /(r) with a modified value of the beam size parameter, 
precisely with vj^ {x) replaced by (uP' (x) /j). Therefore, drS^\x,t) corresponds 
to the expression ( p| ) of dn{x,t) evaluated for a modified beam size: 



dM{x,t)^ 



4 (a2 + ^2^2 



(x) 



exp 



4(^2+^2*2) 



3 J 



(55) 



These expressions are greatly simplified in the limiting case of a small beam 
size {w <^ ar), where dn^^\x,t) is proportional to [x), so that: 



3 



(56) 
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Assuming furthermore that the Rayleigh divergence length is large (cr ^ Ir), 
one may disregard the x-dependence of the on-axis saturation parameter Sm 
and obtain a closed analytical expression for dus {x, t) and therefore for ds'. 

^.W = ^ln(l + 2s„0 (57) 

This expression may then be used to derive the bistability relation which con- 
nects the mean field sent into the cavity and the mean field inside the cavity (see 
for example |9|). An important consequence of this equation is that as {t) and 
(7 (t) have exactly the same time dependence. It follows that the evaluations of 
the number of atoms obtained from the linear and non linear phase shift agree 
as soon as the transverse structure of the mode has been properly accounted 
for. It is however worth to stress that this property relies on the two simplifying 
assumptions w <C Ur and Ur ^ Ir which have been used to derive equation 
(p7[), and does no longer hold in a more general situation. 



5 Fluctuations of the number of atoms 

Up to now, we have studied only the mean value of the field modification in the 
atomic medium. We now want to evaluate the fluctuations of this modification 
which are associated with fluctuations of the number of atoms coupled to the 
field. As shown by the previous discussions, the definition of such a number is 
not simple in the general case. This is why we will study the fiuctuations of the 
number of atoms with the help of simplifying assumptions. 

First, we will restrict our attention to the regime of a linear polarisability 
where the field modification is related to the quantity cr [t] given by equation 
(p3|). We are indeed interested mainly in fluctuations of the linear phase shift, 
which is usually much larger than the nonlinear one. We will then consider that 
the Rayleigh divergence length is much larger than the cloud size {(Jr <C Ir), so 



that the beam area S is independent of x and equation (23) may be rewritten 



in terms of an effective number n[t) of atoms in the probe beam at time t: 



^(t)^^ (58) 



with: 



n{t)= J dn {x, t) = y / (r) p (r, t) dr (59) 

It would be possible to study the fluctuations of a (<) in the general case, using 
the same techniques as presented below. However, it will be more instructive to 
discuss the fluctuations of n (t) , which is a dimensionless number and is expected 
to have a nearly Poissonian statistics. For completeness, note that in the limiting 
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case studied in the present section (cTr <C Ir), n (t) is simply the product of a {t) 
by the beam section S. 

We come now to the discussion of atomic number fluctuations. To this aim, 
we notice that the atomic density p (r, t) may be considered as an average value 
taken over the random distribution of the atomic positions, precisely over the 
random values of the initial positions and velocities: 

p{r,t) = {P{r,t)) (60) 



P (r, i) = ^ <5 (r - r, (i)) = ^ 5 ( r - r, (0) - v, (0) t - -gf (61) 

1=1 i=l ^ ^ 

where i labels the atoms {i = 1 ... A/" with a total number of atoms N in the 
cloud); ri (t) is the position of the ith atom at time t; Vi (0) and v, (0) are its 
initial position and velocity. Therefore, n (t) may be considered as the mean 
value of a random number N (t) : 

Ar(t) = ^^/(r,(i)) (62) 

i=l 

This random number has fluctuations which, as already discussed in the intro- 
duction, may reveal themselves as fluctuations of the linear phase shift and then 
of the cavity-detuning parameter. 

In order to discuss the statistical properties of these numbers, we will con- 
sider that the atoms arc statistically uncorrclatcd in the initial distribution. We 
will define covariances of random variables according to the general prescription: 

{u, v) = {uv) - {u) {v) (63) 

dA/o is the random variable associated with the number of atoms in an elemen- 
tary volume dr dv of the phase space at time t = 0. We consider it to be a 
Poisson random variable: 

(dVo,rfAAo) = (dA/'o) =7r(r,v,0)a'r dv (64) 

while different dMo and dM^ evaluated for non overlapping elementary volumes 
are uncorrelated random variables: 

(dAAo,cW^)=0 (65) 

Before coming to the discussion of the fluctuations of the effective atomic 
number in the detection beam, we want to recall briefly how these assumptions 
are commonly used to derive Poisson statistics. Usually, one studies a number 
of atoms integrated in a finite volume: 
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No - / dNof. (r) 



(66) 



where the function is either 1 inside the detection volume or otherwise; 
this function therefore obeys the property — fv By summing up the vari- 
ances corresponding to the various elementary volumes, one deduces that the 
integrated number Nq has a Poissonian variance: 



{N,,No)= / {dAfo) f^r) = {No) 



(67) 



The same conclusion is reached for the number of atoms integrated in a finite 
volume after a free fall during a time t: 



N{t)^ J dMo fv (^r + vt+ ^gf^ 



(68) 



{N (t) , N (t)) = j (dMo) fl [v + ^t + -gt^ j = {N it)) (69) 

The case studied in the present paper does not lead to such a simple con- 
clusion since the effective number is defined with a Gaussian weight function 
which does not obey the property — f. By summing up the variances cor- 
responding to the various elementary volumes, one obtains the variance of this 
effective number to be: 



{N (t) , N (t)) - j f[v + ^t+ igt^^ IT (r, V, 0) dv dv 



(70) 



By changing variables from {r (0) , v (0)} to {r (i) , v (t)} and integrating over 
the velocity distribution, this variance can be rewritten as: 



(Nit), Nit)) 
while the mean value is: 



fir) p(r,t) dr 



(71) 



(Nit)) 



/(r) p(r,t) dr 



previously for a it) with the effective beam section (5 



TTW 



72): 



(Nit)) 



Ml 



■ exp 



ir^ + rl+t^) 



(72) 



The mean value is obtained as the product of the expression (E3) evaluated 



(73) 



The variance has the same expression evaluated for a modified beam size pa- 
rameter {w"^ replaced by (w;^/2), i.e. replaced by (r^/2)): 
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{N{t),N{t)) 



Ml 



2t2 



2t2 



exp 



2*4 



r| (2r2 



(74) 



' g \" ' r > ' ui ^" 2<2) ^ 

The variance is therefore smaher than the mean value, which implies that the 
distribution of the random variable N (t) has always a sub-Poissonian character: 



{N{t),N{t))<{N{t)) 



(75) 



In the limiting case where the detection beam waist w is much smaller than the 
trap radius ar, the variance is simply half the mean value: 



{N{t),N{t)) = ^{N{t)) 



(76) 



We now look for a dynamical characterization of temporal correlations rather 
than for a static statistical description. Precisely, we wish to evaluate the covari- 
ance function for the numbers of atoms N{t) and N{t/) present in the detection 
beam at two different times. By summing up the variances corresponding to the 
various elementary volumes, one obtains this covariance function as: 



{N (t) , N {t')) = j f[v + wt+ ig<2^ / [t + vt' + IT (r, v,0) dv dv 

(77) 

Since the positions at times t and t' are linear functions of the initial position and 
initial velocity, this expression appears as a six-dimensional Gaussian integral 
which may be readily evaluated. We will write the covariance function: 

C'nn (T, t) = {N (t) , N {t')) (78) 

as a function of the mean time T and of the delay r between the two time 
parameters t and t': 

T='Af r = t-t' (79) 

The time T will be associated with the global time variation of the atomic 
number during the free fall of the cloud while the delay t will rather correspond 
to the correlation between numbers evaluated at different times. 

The explicit evaluation of the Gaussian integral (|7^) gives the following 
expression: 

Cnn (T, r) = n (0) L (T, r) exp [-Af (T, r)] (80) 

where: 

"(0) = r^ (81) 
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L(T,r) = 



M{T,t) = 



(r2 + 




, 2rl) + 4 (r,2 + 4) TV^ 


g 


2r2r2 + 







(82) 



(83) 



The expression of the variance of the number N{T) evaluated at a given time 
T is recovered at t — Q: 



Cnn (T,0) = (iV(r),iV(T)) 



(84) 

Instead evaluating the exact expressions ( p0[]8^ ), we will rather concentrate 
on the limiting case where the detection waist is much smaller than the trap 
radius. This assumption implies in particular that the correlation time is 
much smaller than the global fall time r^. This will allow us to treat the num- 
ber of atoms as a quasistationary random variable. Considering only not too 
long delays r <C and not too short fall times T ^ t^^, we may rewrite the 
covariance function in terms of two scaled time parameters T/tj. and t/t^: 



Cnn (T, T)^n (0) L (T, r) exp hC ("T - brL (T, r))] 



with: 



n(0) 
L{T,r) = 



AT = I — 



5t = 2 ( - 

. Tr 



4+|I 

Tr 



Tr 



^ _9 



(85) 

(86) 
(87) 

(88) 

(89) 

(90) 
(91) 



In the limiting case C 1 where the time scale characteristic of gravity is long, 
the covariance function is essentially a Lorentzian function. More generally when 
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including the effect of gravity, it may be expanded in terms of powers of this 
Lorentzian function: 



°° Tl+k 

Cr,N{T.T)=n{Q)e-^'^-Y.^(:bTt^ (92) 

fe=0 

Since the correlation time, of the order of r^,, is much smaller than the 
characteristic time of variation of the mean values, of the order of t^, the num- 
ber of detected atoms may be considered as a quasistationary random variable 
and it is therefore possible to characterize its fluctuations by a noise spectrum 
Snn iT,uj): 

5Arjv(T,^) = JdT e^'^^C^w(r,r) (93) 

This spectrum is obtained through a Fourier transformation with respect to the 
delay r and it slowly depends upon the time of fall T. It has the same defi- 
nition as the ambiguity function of radar theory [lO| which is the analog for 
time-frequency distributions of the Wigner distribution for position-momentum 
distributions [ pl| . Because of the quasistationary character of the random vari- 
able N (t), it will remain positive for all values of the parameters, as the noise 
spectrum for a stationary random variable. The spectrum may be written as 
the product of the variance at time T by a normalized spectrum: 

Snn (T, lo) = Cnn (T, 0) S (T, to) = (T) Snn (T, lo) (94) 

PSnn{T,uj) = 1 (95) 

In the limiting case ( 1 , the effect of gravity upon the correlation function 
Cnn {T, t) may be disregarded, so that the latter is a Lorentzian function of t 
and the noise spectrum has an exponential shape: 

Snn (T, to) ^ n (0) ^e-"^-""" (96) 
ax 

The normalized spectrum: 

Snn (T, cj) = ^aTT,„e-"^'""- (97) 

has a linewidth {axTw) and a peak value ( ) mainly determined by the 

time of flight through the probe beam. 

In the general case of an arbitrary the noise spectrum may be deduced 
from expression (^) of the correlation function (using formulas 8.432.5 and 
8.468 in ref. [|l|) as: 
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"T t'o V 4 y {k\f 

where pk is a polynomial function of frequency: 
This leads to a normalized spectrum: 



SNNiT,oj) = naTT^e ^ ™e ^ 2^ -p 777-2 (100) 

with a peak value at zero frequency: 

S^^ (T,0) = ^TT.e"? f ^ (101) 

fe=oV 4 y (fc!) 

The expressions ( p8[jl0l| ) provide a quantitative evaluation of the fluctua- 
tions of the number of atoms coupled to the probe beam in the general case 
of a motion under gravity. They may be used for assessing the effect of these 
number fluctuations in any high precision measurements such as spectroscopic, 
nonlinear or quantum optics experiments where the falling cloud acts as the 
atomic medium. For only qualitative assessments, we can draw the following 
two conclusions: first, the instantaneous variance of the effective atomic num- 
ber has a subPoissonian character; it is precisely half the Poissonian variance 
at the limit of a small probe beam size. Second, the correlation time for the 
fluctuations of the effective atomic number is mainly determined by the time of 
flight Tm through the probe beam. 



6 Fluctuations of the cavity detuning 

We now focus our attention on experiments where the probe field goes through 
an optical cavity containing the cloud of cold atoms. As already discussed in 
the introduction, atomic number fluctuations emulate fluctuations of the cavity 
detuning. The present section is devoted to a precise evaluation of this effect. 
For simplicity, we consider only the regime of a linear polarisability in the limit 
where the Rayleigh divergence length is much larger than the cloud size (same 
simplifying assumptions as in section 5). We also consider the dispersive regime 
where the atom is excited far from resonance ((5^ ^ 1). 
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A key parameter in the experiments with atoms inside a cavity is the so- 
.ed cooperativity parameter (see for example [|l3, 14 and references therein): 



called coo_ 

C = (102) 

where k is the decay rate of the intracavity field and Tc the cavity round-trip 
time {2k,Tc is just the intensity transmission of the coupling mirror), while S and 
n are respectively the beam section and the effective atomic number studied in 
the present paper. The effect of the atomic medium on the field may be written 
as a variation $ of the cavity detuning (see for example |p^): 

f . . 14 " (103) 

The fiuctuations of the atomic number N studied in the previous section are 
therefore equivalent to fiuctuations of the cavity detuning $ characterized by a 
noise spectrum: 



that is (using equations ^ and 102): 



C{T) 3A^ SNNiT,u) 
S^^ (T, u) = ^ (105) 

We have introduced here a time-dependent cooperativity parameter C'(T) pro- 
portional to the mean effective number n(T) at time T. Since the cavity detuning 
$ is a frequency, its noise spectrum S^^ als o ha s the dimension of a frequency. 
The first factor appearing in its expression ( |105[ ) is the frequency k , which can 
be regarded as measuring the noise spectrum of detuning fluctuations associ- 
ated with the decay of the intracavity field through the coupling mirror. The 
other factors which appear in equation ( |105| ) have been written in dimensionless 
terms. 

As long as S'$* remains smaller than k, the effect of the detuning fluctu- 
ations on cavity dynamics remains in a linear regime. In the opposite case, it 
would be necessary to give a more elaborate treatment of these fluctuations ac- 
counting for multiplicative noise processing. It is worth stressing that detuning 
fiuctuations may not be considered as a white noise, since the spectrum has a 
width essentially given by the inverse of the time of flight through the detection 
beam. 

In the case where the measurements are performed at frequencies much 
higher than this spectral width, the detuning fiuctuations become negligible. 
This is in particular the case in the experiment on squeezing with cold atoms 

^ where the squeezing spectrum is monitored at frequencies in the MHz 
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range. It follows that the effective atomic number may be considered as a fluc- 
tuationless classical variable, even when it varies as a function of time while the 
cloud is exploding and falling down. The fluctuations of the effective number 
might be detectable at much lower frequencies in the kHz range. In order to 
detect these fluctuations, it would thus be necessary to get rid of the large excess 
noise present in this frequency range. 

In conclusion, we have investigated the problem of evaluating the number 
of atoms effectively interacting with a probe beam passing through a cloud of 
cold atoms released from a magneto-optical trap, accounting for the transverse 
Gaussian profile of the probe beam, in the general case of a motion under gravity. 

We have shown that this problem requires a detailed treatment which de- 
pends on the quantity of interest in the specific measurement. We have given 
explicit evaluations of the effective atomic number for the two cases where the 
linear or nonlinear phaseshifts are measured. These two cases correspond to a 
same time variation of the effective number when the two simplifying assump- 
tions of a small beam waist (w <C dr) and of a long Rayleigh divergence length 
(cTr <C Ir) are fulfilled. The more general situation may also be dealt with by 
using the expressions obtained in sections 3 and 4. 

We have also calculated the correlation functions characterizing the fluc- 
tuations of the effective atomic number, focussing attention upon fluctuations 
of the linear phaseshift since they may be quite large and consequently affect 
high precision measurements such as spectroscopic, nonlinear or quantum op- 
tics experiments. We have found that the instantaneous variance of the effective 
atomic number has a subPoissonian character, due to the Gaussian profile of 
the probe beam. In particular, the variance is half the Poissonian variance at 
the limit of a small beam waist (w ^ tXr). At the same limit, the spectral width 
characteristic of these fiuctuations is essentially the inverse of the time of flight 
through the probe beam for an atom flying at the mean thermal velocity. The 
effect of these number fluctuations on the field fiuctuations is thus confined to 
a low frequency range, especially for cold atoms which correspond to relatively 
long time of flight. For atoms at room temperature, the spectral width would be 
larger. In both cases, the techniques used in the present paper allow to obtain 
quantitative estimates of the influence of atomic niunber fluctuations on field 
fluctuations, and then to the sensitivity of optical measurements. 
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